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Abstract 

The theory of the 'proton spin' effect proposed in our earlier papers is extended to 
include the chiral SU(3) symmetry breaking and flavour mixing induced by non- vanishing 
quark masses in QCD. The theoretical basis is the derivation of exact, unified Goldberger- 
Treiman (GT) relations valid beyond the chiral limit. The observed suppression in the 
flavour singlet axial charge a°(Q 2 ) is explained by an anomalously small value for the 
slope of the singlet current correlation function (0\T d^J^ 5 <9^J° 5 |0), a consequence of the 
screening of topological charge in the QCD vacuum. Numerical predictions are obtained 
by evaluating the current correlation functions using QCD spectral sum rules. The results, 
a°(Q 2 ) = 0.31 ± 0.02 and / dx g[{x, Q 2 ) = 0.141 ± 0.005 (at Q 2 = 10 GeV 2 ), are in good 
agreement with current experimental data on the polarised proton structure function g\. 
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1. Introduction 

The 'proton spin' problem, i.e. the question of why the first moment of the flavour 
singlet component of the polarised proton structure function g\ is anomalously suppressed, 
has inspired an impressive research effort, both theoretical and experimental, for over a 
decade. (For recent reviews, see e.g. refs.[l,2]. As is well-known, the first moment of g\ 
can be expressed in terms of the axial charges of the proton as follows: 

j\xg{{x-Q 2 ) = lcf s (a s (Q 2 )) ( a 3 + V) + ±C7? (a s (Q 2 )) a°(Q 2 ) (1.1) 

Here, Cf s , Cf are the appropriate Wilson coefficients arising from the OPE for two 
electromagnetic currents, while the axial charges are defined from the forward matrix 
elements 

(p,s\J%\p,s) = \a 3 ^ (p,s\J 8 5 \p,s) = (p, s\ J%\p, s) = a°(Q 2 )s M (1.2) 

where J^ 5 (often denoted Afy are the axial currents and s^ is the proton polarisation 
vector, a 3 and a 8 are known in terms of the F and D coefficients from beta and hyperon 
decay, so that an experimental determination of the first moment of g\ in polarised deep 
inelastic scattering (DIS) allows a determination of the singlet axial charge a°(Q 2 ). The 
'proton spin' problem is that it is found experimentally that a°(Q 2 ) is strongly suppressed 
relative to a 8 , which would be its expected value if the OZI (Zweig) rule were exact in this 
channel. 

DIS is normally described theoretically using the QCD parton model. In this model, 
the axial charges are represented (in the AB renormalisation scheme) in terms of moments 
of parton distributions as follows[3,4]: 

a 3 = Au-Ad a 8 = Au + Ad-2As a°(Q 2 ) = Au + Ad+ As-n f ^Ag(Q 2 ) (1.3) 

In the parton model, the 'proton spin' problem takes the following form. In the naive, or 
valence quark, parton model we would expect the strange quark and gluon distributions 
to vanish, i.e. As = 0, Ag(Q 2 ) = 0. In that case, a = a 8 , the OZI prediction. Inserted 
into eq.(l.l), this gives the Ellis-Jaffe sum rule [5]. However, the observed suppression 
a°(Q 2 ) < a 8 can be accommodated in the full QCD parton model by invoking either 
or both a non-zero polarised strange quark distribution As ^ or a non-zero polarised 
gluon distribution Ag(Q 2 ) ^ 0. An interesting conjecture (in line with the insights of our 
alternative approach) is that the suppression is primarily due to the gluon distribution[3], 
although a quantitative prediction would still only follow if Ag(Q 2 ) can be independently 
measured, either through a precise analysis of the Q 2 dependence of g\ [4] or directly 
through other less inclusive high energy processes such as open charm production. Notice, 
however, that even in the QCD parton model picture, it is not possible to identify a (<5 2 ) 
with spin[6]. This identification only holds for free quarks, in which case the Q 2 scale 
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dependence (which is related to the {7a (1) axial anomaly) disappears from a . As has 
been emphasised many times[7-9,l], the so-called 'proton spin' problem is not a problem 
of spin - rather, it is a question of understanding the dynamical origin of the OZI violation 

a°(Q 2 ) < a 8 . 

In a series of papers[7, 10-12], we have a proposed an alternative, complementary, 
approach to the 'proton spin' problem which provides both a new physical insight into 
the origin of the suppression in a°(<5 2 ) and a quantitative prediction, which is in good 
agreement with the current experimental data. In our approach (reviewed in refs.[8,l]), 
the flavour singlet axial charge a°(Q 2 ) is decomposed into the product of an RG-invariant 
1PI vertex and a non-perturbative, but target independent, QCD correlation function 
which we subsequently identify as the first moment of the QCD topological susceptibility 
X '(0),i.e. 

a°(Q 2 ) = r^ev^f^ (1.4) 
where x'(0) = ^x{k 2 )\ k=0 and 

X (k 2 ) =ijd 4 x e lkx (0\T Q{x) Q(0)|0> (1.5) 

with Q = f^trG ^vG^ the gluon topological charge density. In the vertex, rf denotes 
the 'OZI Goldstone boson', i.e. the (unphysical) state which would become the Goldstone 
boson for spontaneously broken £/a(1) in the absence of the anomaly (OZI limit). We now 
make the key assumption that the RG-invariant vertex obeys the OZI rule, viz. T v »nn = 
s 7v iv- Then, comparing with the GT relation for a 8 in the chiral limit, we obtain our 
prediction: 

T - { 

The suppression in a°(Q 2 ) simply reflects an anomalously small value of x'(O), which we 
confirmed by an explicit calculation using QCD spectral sum rules[12]. In this picture, 
the suppression is therefore a target-independent effect, i.e. a generic property of the QCD 
vacuum rather than a specific property of the proton. Since a°(Q 2 ) can be expressed, via 
the axial {7a (1) anomaly, as the proton matrix element of the topological charge Q, its 
suppression can be understood as a screening of topological charge in the QCD vacuum. 
We conclude that the fundamental dynamics underlying the violation of the Ellis-Jaffe 
sum rule is not to do with quark spin, but instead is a manifestation of topological charge 
screening by the QCD vacuum. 

This 'target-independence' property may in principle be tested directly in experiment, 
by studying semi-inclusive DIS in which a hadron carrying a large fraction of the target 
energy is observed in the target fragmentation region. Such experiments should be possible 
at, e.g. polarised HERA. The details of our proposal are described in refs. [13-15]. 

Although our previous analysis was restricted to the chiral limit of QCD, it was clear 
from the derivation that the final result for a°(<5 2 ) should only have a weak dependence 
on the quark masses. Indeed, the cancellation of the explicit quark mass dependence in 
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a°(Q 2 ) was already pointed out in ref.[7]. However, the introduction of quark masses does 
involve a significant complication in the analysis. The purpose of this paper is therefore 
to present the generalisation of our theory of the 'proton spin' to QCD beyond the chiral 
limit. 

Underlying our approach to polarised DIS - the 'CPV method reviewed in some detail 
in ref.[l] - is the C/a(1) Goldberger-Treiman (GT) relation. The initial insight that the 
'proton spin' problem was essentially one of OZI violation and could be understood in 
terms of an extension of the GT relation to the flavour singlet channel was made by one 
of us in ref.[7]. The Ua(X) GT relation was subsequently put on a firm field-theoretical 
basis in refs. [10,11], where the connection with the slope of the topological susceptibility 
was first established. 

The introduction of quark masses complicates the £7a(1) GT analysis in two ways. 
First, the pure gluon topological susceptibility has to be combined with correlation func- 
tions of pseudoscalar operators of the form 4>$ = <?75<? arising from the extra explicit 
chiral symmetry breaking terms in the divergence of the axial current. Second, since the 
quark masses break flavour SU(3), the UaO-) GT relation is mixed with those for the 
flavour non-singlet axial charges. In this paper, we present the unified GT relations for 
the flavour singlet and non-singlet axial charges, expressing them in terms of 1PI vertices 
and the first moments of the correlation functions of the divergences of the axial currents 
as follows: 

G°A = 7t Fab ^ri b NN (1-7) 

where a = 0, 3, 8 is the SU (3) flavour index, the axial charges are normalised as 

G\ = \a* G a = ^ G A (Q 2 ) = a°(Q 2 ) (1.8) 

and 

(FF T ) ab = hm ^ ij d * x (0 |T d^J^(x) <9^ 5 (0)|0> (1.9) 

It is important to realise that, like the singlet Ua(X) GT relation, these are exact relations 
in QCD field theory. They go beyond the pole-dominance approximation used in the 
conventional GT relations, replacing the decay constants by current correlation functions 
and the ttNN or r/NN couplings by appropriate 1PI vertex functions. 

In the second part of the paper, we make this analysis quantitative by calculating 
the relevant current correlation functions using QCD spectral sum rules. This confirms 
the stability of our earlier results [12] to the introduction of quark masses. A very care- 
ful analysis is presented, with particular attention given to the inclusion of higher order 
corrections in a s or m 2 s r and to the stability of the sum rules with respect to the Borel 
parameter r. We find that the size of the explicit chiral and SU(3) breaking effects is in 
line with our theoretical expectations and confirms the validity of the Laplace sum rule 
method for calculating correlation functions, such as the topological susceptibility, in the 
flavour singlet channel. 

The paper is organised as follows. In section 2, we quote the (anomalous) chiral 
Ward identities which we need to derive the unified Goldberger-Treiman relations. The 
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derivation itself is presented in section 3. In section 4, we relate these new GT relations 
to the 'proton spin' sum rule and describe our expectations for the results of a spectral 
sum rule calculation of the relevant current correlation functions. This is elaborated on in 
section 5, where the pattern of explicit mass cancellations is exhibited in the framework of 
an effective lagrangian. 

In section 6, we evaluate the current correlation functions (1.9) using the QCD spectral 
sum rule technique. The presence of quark masses requires a substantially more compli- 
cated analysis than we had previously performed in the chiral limit to evaluate the slope 
of the topological susceptibility. However, the results confirm our earlier conclusions [12]. 
In section 7, these numerical results are used in the GT relations and our prediction for 
the 'proton spin' suppression is obtained. Our conclusions from this work are presented in 
section 8. 

Appendix A deals with some technical aspects of the derivation of the unified GT 
formulae. Their relationship with conventional current algebra and PCAC methods is dis- 
cussed briefly in Appendix B. The renormalisation group properties of the Green functions 
and vertices involved in the GT relations are derived in Appendix C. Finally, since both 
our formal methods and spectral sum rule analysis have been repeatedly criticised in the 
literature by Ioffe (see e.g. refs.[16,17] and references therein), we explain in Appendix D 
why these criticisms are not correct and show in detail the errors in refs. [16-20] which lead 
Ioffe to his false conclusions. 



2. Chiral Ward Identities 

The derivation of the generalised Goldberger-Treiman relations is based on the chiral 
Ward identities satisfied by the composite operator propagators and vertex functions. In 
this section, we derive these in the form which will be most convenient for the applications 
which follow. 

For the propagators, the starting point is the Ward identity for the generating func- 
tional VF[V^ 5 , V®, 6>, S$, S a ] of Green functions which are '1PI' with respect to the desig- 
nated fields (composite operators). Here, V® 5 , V®, 9, S§ , S a are the sources for the com- 
posite operators J^ 5 , J^, Q, (fi a respectively, where 

Jfc = Tivl$T a q J« = q ltl T a q Q = ^trC^C" 

cf> a 5 = q l5 T a q <p a = qT a q (2.1) 

In this notation, T l = are flavour SU(nf) generators, and we include the singlet £/a(1) 
generator T° = 1 and let the index a = 0,i. We will only need to consider fields where i 
corresponds to a generator in the Cartan sub-algebra, so that a = 0, 3, 8 for rif = 3 quark 
flavours. We define <i-symbols by {T a , T b } = d a b c T c . Since this includes the flavour singlet 
U(l) generator, they are only symmetric on the first two indices. For rif = 3, the explicit 
values are cfooo = ^033 = ^088 = 2,rf 3 3 = d 8 so = 1/3,^338 = <^383 = -^888 = l/v 7 ^. (For 
further notation and description of the formalism used here, see ref.fll].) 
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The chiral Ward identities are written in the functional formalism as follows: 



SW n c SW , d 5W , ,<W , ,<yw 

" 2w /*«o-^- - ^ dc m d — + cU^ — - dadcSl— = (2.2) 



where we have displayed both the anomalous breaking term for the £/a(1) current and the 
soft breaking induced by the quark masses. The quark mass matrix is written as m a T a , 
so that for to/ = 3, 



m u 





° \ 







m d 


=m°l+m s T 3 + m 8 T 8 


(2.3) 








m s J 





Notice that there are no variation terms for the currents themselves as we have restricted 
to fields where a is a Cartan sub-algebra index. It is convenient to use a notation where a 
functional derivative is denoted simply by a subscript. So, also transforming to momentum 
space, we can rewrite eq.(2.2) compactly as 

ifc M W v# - - 2n f 6 a0 W e - d adc m d Wsi + d adc S d W st - d adc S%W S c = (2.4) 

The Ward identities for composite operator Green functions are derived by taking 
functional derivatives of this basic identity. We will need the following identities for 2- 
point functions: 

ikuW V a V b -2rif8 a oW 9V b - M ac W S c V b =0 
ikfjWy^e - 2nf5 a0 W ee - M ac W s -e = 

ik^Wya gb - 2n f S a0 W esb - M ac W s?S b - $ ab = (2.5) 
where we have introduced the still more compact notation 

M ab = d acb m c $ a6 = d abc {(j) c ) (2.6) 

(4> c ) is the VEV (qT c q), so that for to/ = 3 the condensates may be written as 



(mm) 
(dd) 
(ss) 

Notice that with these definitions, 



0' (dd) ] = ^°)l + 2((/) 3 )T 3 + 2((/) 8 )T 8 (2.7) 



-detM = m u m d m s -det<E> = (uu) (dd) (ss) (2.8) 

8 6 

with the obvious generalisation for arbitrary to/. 

Combining the individual equations in (2.5), we find the important identity: 

k^Wya^ - M ac $ cb = W (2.9) 



where Sp is the source for the current divergence operator D a = 2nf5 a oQ + d ac dm c (j)^. In 
conventional notation, 

W S a SbD = tjdx e ik * (0|T D a (x) D b (0)\0) 

= i J dx e lkx (0|T d^J^x) d u 4 5 {0)\0) (2.10) 

The zero-momentum Ward identities play a special role. These follow immediately 
from eqs.(2.5) under the assumption that there are no massless particles (in particular, 
no exact Goldstone bosons) contributing 1/k 2 poles in the 2-point functions. With this 
assumption, we find simply 

2n f S a0 Wee + M ac W Si e = 

2n f 8 a0 W eS b + M ac W S o S b + $ a6 = (2.11) 

The generating functional for proper vertices, T, is defined from W by a partial Leg- 
endre transform (Zumino transform[21]), in which the transform is made only on the fields 
Q, </>5, 4> a and not on the currents. A number of important results on these transforms are 
collected in Appendix A. The resulting proper vertices are '1PI' wrt the propagators for 
these composite operators only. As explained fully in ref.fll], by separating off the particle 
poles in the propagators, this is the definition which gives the closest identification of these 
field-theoretic vertices with physical low-energy couplings such as g n NN etc. 

We therefore define the generating functional r[V^ 5 , V£, Q, 4> a ] as: 

r[y; 5 , y;, q, c <p a ] = w[v^ v?, e, s%, s a ] - J dx (oq + s^ a 5 + s a r) (2.12) 

The chiral Ward identities corresponding to eq.(2.4) are therefore: 

ife M r v - H - 2n f S a0 Q - d acd m c (j)i + d acd (j) d r n - d acd $jl> = (2.13) 

The Ward identities for the 2-point vertices will also be important in the derivation of the 
GT relations. These follow directly from eq.(2.13): 

iknTya yb + a cF frcyb = 

ik^Ty^Q - 2n f 5 a0 + $ ac r^g = 

ik^Tya^b + ^acV^b -M ab = (2.14) 

It is then straightforward to derive the following important identity, analogous to eq.(2.9): 

k^k v Tya yb + M ac $ cb = $ ac r\c^$ d6 (2.15) 
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3. The Goldberger-Treiman Relations 



We now present a unified derivation of the Goldberger-Treiman relations for the flavour 
singlet and non-singlet axial charges. ( 1 ) The derivation follows the principles of refs. [10,11], 
extended to include non-zero quark masses and flavour mixing. This allows us to present 
all the GT relations in terms of a single, unified formula involving the 2-point correlation 
functions of the divergences of the flavour singlet and non-singlet currents. 

The axial charges G\ are defined as the form factors in the forward nucleon matrix 
elements of the axial currents, viz. 



(p, s\J^\p, s) = G\s^ 



(3.1) 



where p^ and s M = u(p, 3)7^75^(25, s) are respectively the momentum and polarisation 
vector of the nucleon. 

To express this matrix element in our composite propagator-vertex formalism, we 
introduce an interpolating field N and source Sn for the nucleon as in refs. [10,11]. (Notice 
that this is purely a formal device - there is no dynamics implicit in this manoeuvre.) The 
matrix element is then just the 3-point function Wy^ B s N s N with the external propagators 
amputated. From eq.(A.8), we can re-express this in terms of the vertex functional V as 
follows: 



(P, s\J^ 5 \p, s) = u(p,s) 



u(p, s) 



W 



SnSn 



VfcSnS* W S N S N 



u{p, s) 



Fv« 5 NN + W V a r j T QNN + Wya^gb T (j} b NN 



u(p,s) (3.2) 



Since the propagators on the rhs vanish at zero momentum (this requires the absence of 
any 1/k 2 poles, which is assured by the Ua(1) anomaly and quark masses), we find simply 



2mNG a A wy$u = u 



u 



(3.3) 



where uln is the nucleon mass. The GT relations then follow immediately from the Ward 
identity (2.13) for V. Differentiating wrt the nucleon fields, we find 



2m N G a A = $ a bT jbxxl^Q 



(3.4) 



where for convenience we define iu T^nn u = T^nn ^75^, etc. 



Similar generalised GT relations have been presented some time ago, with the same motivation of 
explaining the 'proton spin' problem, by Efremov, Soffer and T6rnqvist[22] using the more conventional 
PCAC language. 
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A non-forward version of the GT relation, which was used extensively in refs. [10-12], 
can be found from eq.(3.2) by using the Ward identities (2.5) for the propagators together 
with (2.13) for Ty^NN- This allows us to write, for all k, 

2m N G a A (k 2 ) + k 2 G a P (k 2 )=-(2n f 5 a oWe e +M ac WsiB)fQNN 

- (2nf6 a0 W es > + M ac W SiS ^f^ NN (3.5) 

Gp(k 2 ) is the pseudoscalar form factor in the non-forward matrix element (p, s\J%\p, s), 
and again has no 1/k 2 pole for the reasons given above. This expression clearly reduces to 
eq.(3.4) on using the zero-momentum Ward identities (2.11) for the propagators. 

The remaining step to convert eq.(3.4) into the useful form of the GT relations is to 
normalise the field 0f appropriately. Clearly, eq.(3.4) is independent of the normalisation. 
However, with a suitable choice, the vertices can be made both RG invariant and essentially 
identical to the physical Goldstone boson couplings g^NN etc. To achieve this, we define 
normalised fields 

V a = Ba^l (3.6) 
where B is a constant matrix^ 2 ) such that 

o = 6 ab (3.7) 

This condition ensures that the fields rf have unit coupling to the Goldstone bosons. The 
case of the singlet rf is of course special, since it is only after mixing with the topological 
field Q (and then flavour mixing) that it becomes the physical rj'. The intricacies of this 
are discussed fully in ref.[ll]. Indeed, this is why it is most convenient to impose the 
normalisation condition as above on the matrix of 2-point vertices r^b, which is the 
inverse of the pseudoscalar propagator matrix, since this most simply characterises the rf 
before mixing with Q. 

( 2 ) From the Ward identity (in matrix notation) 

ik^r v ^ v + <t>B T r vv - MB- 1 = 

and writing 
we find 

T vv = k 2 (^B T )- 1 f + (<PB T )- 1 M<S>(B<S>)- 1 

The normalisation condition (3.7) can therefore be applied while keeping B in eq.(3.6) a constant. This 
determines (see Appendix B for the relation to PCAC) 

/(0) = $B T 

However, it is not possible, as in the chiral limit case[ff], to impose a normalisation valid for all k, such 
as r^a^b = k 2 5 a b — (m 2 ) a b, by allowing B to be a function of k 2 . 



—r 
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The proof that the vertices T v a NN defined with the fields normalised according to 
eq.(3.7) are RG invariant now goes through in the same way as shown in the chiral limit 
in ref.[ll]. In Appendix C, we summarise some of the main results for the extension to 
non-zero quark masses. Ref.[ll] also contains a careful description of how our formalism 
is related to standard current algebra (PC AC). A simple illustration of this is included in 
Appendix B. 

Re-expressing eq.(3.4) in terms of the properly normalised vertices, we have 

2m N G a A = $ ac S4 

1 r] b NN | 

(3.8) 

where B is to be determined from the 2-point vertex condition 

^2 F <^L =0 = B ac jfr2 T V c V d \k=0 Bdb = B ac B cb (3.9) 

The straightforward approach to finding T^a^b is to write it as one component of the 
inverse of the propagator matrix Wss (with S = {6>, S§}). This was the approach used in 
refs. [10,11] to relate T^o^o to the topological susceptibility Wee- However, inverting this 
matrix in the multi-flavour case is cumbersome, so here we use an alternative approach. 

First, we use the Ward identities (2.14) to express T^^b in terms of T v ^ V b^: 

kuk u Tya yb = ik u Q ac T JyCyb 
At5 c5 ^5 vo 



= $ ac V^a $ db - M ac $ cb (3.10) 

Then, using the general identity (A. 9) for partial Legendre transforms, we relate Ty^y^ 
to the 2-current propagator W V a V b : 

Tya yb = Wya yb — Wya 5 W7 T Wn-yb (3-H) 

where S and T represent {6>, S^}. Combining with eqs.(2.9) and (2.15) finally yields 

$ac r^c^d $ db = W S a DS b D + k^Wya^s W^r W TV bk v (3.12) 
and so, in matrix notation, 

® BTm = ^ W s D s D \ k=G + ^(M^ 5 s W sr Wrv v5 K)\ k=0 (3.13) 

The argument is almost complete. $S T is precisely the combination we need for the 
GT relations (3.8), and is related by eq.(3.13) directly to the first moment of the 2-point 
correlation function (2.10) for the divergences of the currents. This generalises the relation 
with the topological susceptibility found in refs. [10,11]. 

It remains to show that the final term in eq.(3.13) vanishes. The first and last fac- 
tors are of 0(/c 2 ), so this will contribute zero unless there is a 1/k 2 pole in the inverse 
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pseudoscalar propagator matrix Wg-^. As already mentioned, there are no 1/k 2 poles in 
the propagators themselves, so all we need show is that the determinant A(k 2 ) of this 
propagator matrix is no n- vanishing at k = 0. This follows from the formula 

A(0) = Wee(0) (detM)- 1 det$ (3.14) 

since detM and det$ are non-zero (eqs.(2.8)) and Wgg is non- vanishing away from the 
chiral limit. An elegant proof of eq.(3.14), based on the zero- momentum Ward identities 
for Wst is as follows. Define 

with A = detW. Then, 

detM A = detM detW 



W, 



0S£ 



M ac W s? e M ac W S o S t 



(3.16) 



Using the zero-momentum Ward identities (2.11) and, in the case of the a = row taking 
a linear combination with the first (6) row, the determinant simplifies, leaving 



detM A = 



-$ab 

= Wee det$ 



(3.17) 



as required. 

This completes the proof of the GT relations. To summarise, we have shown that the 
flavour singlet and non-singlet axial charges are given by the single, unified relation 

2m N G a A = F ab r v b NN \ k=Q (3.18) 
where F = $S T is determined from 

FFT = ^2 W s»s D \ k=0 (3-19) 

that is, 

F ac Fj b = \im^ijdx e ikx (0|T d»J^(x) a"J* 5 (0)|0> (3.20) 

In current algebra terms, the matrix F determined from eq.(3.20) can be identified, 
subject to the standard PCAC (pole dominance) approximation described in Appendix B, 
with the pseudo-Goldstone boson decay constants. Under the same PCAC assumptions, 
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the vertices can be identified with the low-energy meson-nucleon coupling constants, i.e. 
T v 3 NN\k=o = 9-kNN-, etc. It should be emphasised, however, that the relations (3.18) and 
(3.20) are exact - they are neither dependent on this interpretation nor make any PCAC 
assumption or approximation. 

Eq.(3.18) is therefore the natural generalisation of the U(l) GT relation proved in 
refs. [10,11]. There we showed that, in the chiral limit, the singlet axial charge is given by 

2m N G° A = 2n f ^/^j f v o NN \ k=0 (3.21) 
where %'(0) is the first moment of the topological susceptibility, viz. 

X '(0) = )im-^ijdx e ikx {0\T Q{x) Q(0)|0> (3.22) 

The new formula extends this to include mixing with the flavour non-singlet sector, intro- 
ducing a matrix structure which replaces the simple square root in eq.(3.21) and general- 
ising the fields in the correlation function to the entire divergence of the current including 
mass terms as well as the anomaly Q. 



3.1 Flavour mixing 

The remaining theoretical question related to the unified GT relation (3.18) concerns 
SU(3) flavour mixing and the extent to which eq.(3.20) determines F. 

In the context of PCAC or chiral perturbation theory (incorporating the 1/N C expan- 
sion to allow the inclusion of the flavour singlet rj' (see ref.[23] and references therein)), 
the problem of 77 — 77' mixing is currently receiving renewed attention[23,24]. According 
to Kaiser and Leutwyler[23], the decay constant matrix in the a = 0, 8, rj — rj' sector is 
required to contain 4 parameters, viz. 

(fori' fo v \ _ /Vocostfo -iosini?o\ / 3 23 n 

V hv' fori J \h sin ^8 /8 COS t?8 J 

with i) 7^ i? 8 , in contrast to the previous conventional analysis which assumed $0 = "#8- 
Low-energy theorems[23] yield /o, fs and sin(??o — ^s) easily, from the diagonal and off- 
diagonal combinations J2p= v ',ri /op/op, J^pfspfsp and J^pfopfsp respectively, but the 
remaining combination of the mixing angles is harder to identify. 

Returning to our result (3.18), (3.20), we find a similar situation. Clearly the relation 
(3.20) can only determine 3 parameters in F a b (in the a,b = 0,8 sector). However, in 
general our analysis requires F a b to be characterised by 4 parameters, including a mixing 
angle undetermined by eq.(3.20). 

To see this in more detail, consider the defining equation (3.9) for B, re-expressed in 
terms of F itself: 

$ « c ^2 r ^^L=0 $dfe = Fac j^ T V c V d \k=0 F db = F acFj b (3.24) 
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Since the l.h.s. is symmetric, it can be diagonalised by an orthogonal matrix R, i.e. 




= R T D 2 R 



(3.25) 



where D is diagonal. But since the r.h.s. of eq.(3.24) is unchanged if F is right-multiplied 
by an independent orthogonal matrix O t , we find the general solution 



The four parameters are interpreted as an overall scale / which in the PCAC approximation 
and the limit of exact SU(3) becomes the pion decay constant, an OZI breaking (s 7^ 1) 
parameter, an angle 9 characterising SU (3) breaking, and finally a mixing angle cf> for the 
rf fields. 

We can now check whether any special cases of this general parametrisation are con- 
sistent with the RGE (C.10) for F a b. Clearly there is an exact SU(3) limit in which both 9 
and 4> are zero, with the RGE being satisfied by Vs = 7s, Vf = 0. The case 9 = 0, 7^ 
is also RG consistent, with the same solution. (This corresponds to the case i} = $8 m 
the parametrisation (3.23).) However, the RGE cannot be satisfied with 9 7^ but = 0. 
In the general 577(3) breaking case, therefore, F a b must depend on a mixing angle <p which 
cannot be determined from the condition (3.20). 

The consequence of this is that in order to use the full unified GT relation (3.18) to 
make predictions that include the effect of SU (3) flavour mixing, we would need a further 
(linear) condition on F beyond eq.(3.20). 

Without such an extra condition, we have to neglect flavour mixing. This will produce 
an uncertainty in our final predictions which we can estimate by defining an SU(3) breaking 
parameter from the condensates (2.7), viz. t = V6(0 8 )/(0°). (The y/6 factors here and 
in the definition of s arise because of the different normalisation of the singlet generator: 
trT°T° = rif whereas trT^T- 7 = for i,j = 1, . . .8.) Using standard values (6.16) for 
the condensates, t ~ 0.16. We therefore expect the 577(3) breaking angle sin# to be of 
0(t) and so omitting flavour mixing effects will produce an uncertainty of O(10 — 20%) in 
our final results. 

However, when we evaluate the correlation functions (3.20) using QCD spectral sum 
rules (section 6), it is in any case difficult to do better, since to incorporate flavour mixing 
we would have to saturate the spectral functions (below t c ) with the two states rj' and 77, 
for both the flavour diagonal and off-diagonal correlators. This would greatly complicate 
the analysis without significantly improving the accuracy of the final results, so we do not 
attempt it here. Consequently, our final predictions, presented in section 7, are subject to 
the approximation of neglecting SU (3) flavour mixing. 



F = R T DO 



(3.26) 



A convenient parametrisation for F is therefore 




(3.27) 
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4. GT Relations and the Sum Rule for the First Moment of g\ 



As already quoted in the introduction, the first moment of the polarised proton struc- 
ture function g\ satisfies the following sum rule: 

Y\{Q 2 )^ C dx g\{x-Q 2 ) 
Jo 

= ^Cf s (a s (Q 2 )) (a 3 + la 8 ) + ^ (a s (Q 2 )) a°(Q 2 ) (4.1) 
where the a a are the axial charges occurring in the GT relations we have just derived: 

G\ = \a* ^ = ^7r 8 G\{Q 2 )=a\Q 2 ) (4.2) 

The RG scale (Q 2 ) dependence of a°(Q 2 ), due to the anomalous dimension of the singlet 
axial current, is explicitly displayed. 

This sum rule has been analysed using the composite operator propagator-vertex 
approach to DIS in ref.[12]. Working in the chiral limit, it was shown that 

TUnglet = \^ 2n f Cf {a s (Q 2 )) V7(0j\ Q 2 f>JViV (4-3) 

Following refs. [10-12], we now assume that the vertices are well approximated by their OZI 
values. This is the key assumption that allows us to make a quantitative prediction for 
on the basis of a calculation of the topological susceptibility alone. The RG invariance 
of the vertices is a necessary condition for this assumption to be reasonable. Further phe- 
nomenological evidence from Ua(1) current algebra supporting this conjecture is discussed 
in refs. [10-12]. 

We therefore assume that T v o NN satisfies the OZI rule, i.e. T v o NN = y/2 T v s N n, 
so that all the OZI breaking in T\ sing i et resides in the topological susceptibility a/x'(0). 
Comparing with the standard OZI relation for a 8 , we then find that [10, 11] 

a°(Q 2 ) V6 



^ = ^ V»| Q2 (4.4) 

This leads to the following prediction[12] for the 'proton spin' sum rule: 

a°(Q 2 = lOGeV 2 ) = 0.35 ± 0.05 

Y\{Q 2 = lOGeV 2 ) = 0.143 ±0.005 (4.5) 
based on our original derivation of x'(0) using QCD spectral sum rules[12]: 



V / X 7 (0)| Q2=10 GeV 2 = ( 23 - 2 ± 2 - 4 ) MeV ( 4 - 6 ) 
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This is to be compared with the OZI prediction a \ OZI = or = 0.58 ± 0.03 and with 
the current experimental data from the SMC collaboration [25]: 

r?(Q 2 = 10GeV 2 )| = f dxg p 1 (x;Q 2 ) = 0.141 ±0.012 (4.7a) 

V ' ' 7o.003 

The result for the entire first moment depends on how the extrapolation to the unmeasured 
small x region x < 0.003 is performed. Using a simple Regge fit, SMC find F p = 0.142 ± 
0.017 from which they deduce a = 0.34 ±0.17, while using a small x fit using perturbative 
QCD evolution of the parton distributions [4,26] they find T p = 0.130 ± 0.017 and a = 
0.22 ± 0.17 (all at Q 2 = lOGeV 2 ). 

More recently, SMC have published a further analysis[27] of their data, this time 
quoting a slightly lower number for the integral over the measured region of x: 

T P (Q 2 = lOGeV 2 ) | (x>Q 0Q3) = J ^ dx g{{x- Q 2 ) = 0.133 ± 0.009 (4.76) 

Whatever the ultimate resolution of the small x extrapolation, it is encouraging that 
our prediction (4.5) is in the region now favoured by the data. We would therefore like 
to test the precision of our result further. An assumption in making this prediction was 
that the value of a is smooth in the quark masses, so that the chiral limit will be a good 
approximation, correct up to the usual order of soft 577(3) breaking in ratios of decay 
constants. The extended version of the GT relations derived in section 3 allow us to test 
this. 

Using the new form (3.18) of the GT relations, we can immediately rewrite the first 
moment sum rule for g\ in the following compact form: 

Y\{Q 2 ) ee J\x g{{x,Q 2 ) = i JL C<?(a s (Q 2 )) G\ 

= Cl{a s {Q 2 )) F ab t vbNN (4.8) 

where we have defined the numerically rescaled Wilson coefficients = Cf and Cf = 
V3Cf = ^C^ s . All the vertices T v a NN are RG scale invariant. Apart from the running 
coupling in the Wilson coefficients, the only Q 2 scale dependence is contained in the singlet 
components of F a &, according to the RG equation (C.10). 

Following the approach of refs. [10-12], we again assume that these RG-invariant ver- 
tices are all well approximated by their OZI values. Isospin invariance ensures that the 
off-diagonal terms in the 2-current correlation functions involving the triplet index 3 vanish, 
so we need consider only: 

d ( {^J% ^J° 5 |0) <O|0" J% 0"J* B |O) \ 

FF T = lim —2 i J (0\d»J* 5 PJ* 6 \0) (4.9) 

k -° dk V<°l^5^7° 5 |o> o (o\d»4 5 d»j* 5 \o)J 
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In section 6, we evaluate these correlators (in the 0, 8 sector) using spectral sum rules. As 
explained in the last section, we restrict ourselves to single-particle saturation of the appro- 
priate spectral functions, keeping only the rj' contribution in the flavour singlet correlator 
and rj in the octet. We do not evaluate the off-diagonal correlator, which is expected to 
be small due to the approximate cancellation between the decay constants for the rj' and 
rj states in the spectral function. 

In this approximation of consistently neglecting the SU(3) flavour mixing, and using 
the OZI relation T v o NN = T v s NNl we then identify the 'proton spin' suppression from 
eq.(4.8) as: 



a 8 v/6 F 88 



where 

Fqo 

and 



lim A z(0|^j0 5 ^J° 5 |0> (4.11) 



? 88 * ^/hm A z(0 |aMj8 5 ^J8 5 | ) (4.12) 



5. Quark Mass Dependence 

Before presenting the spectral sum rule analysis in section 6, we give here some ana- 
lytic arguments based on the Ward identities which suggest that the first moment of the 
correlation functions (4.9) will have only a weak dependence on the quark masses. The 
basic reason for this is the identification, up to the standard PCAC pole-dominance as- 
sumptions, of F with a decay constant. The explicit dependence on the pseudo-Goldstone 
boson masses which is present in Wga^ drops out in the first moment. 

To see this another way, consider the Ward identity (2.9). Taking the first moment, 
we see that the explicit dependence on the quark masses m vanishes, leaving 



FF T = -!Lwaao b \, „ = -?L(k»k v W va 



dk 2 ^^Ito dk 2 V V^s 

The 2-current correlation function can be parametrised as 



k=0 



(5.1) 



WV- B v«- 8 = ijdx e ik ' x (0|T Jfe(x) J* 5 (0)|0) 



= n^ 2 )(^-^) + uf(k 2 )^ (5.2) 

where IIj^/c 2 ), n^/c 2 ) are dynamical functions which are not determined by the chiral 
Ward identities. (They are of course related by eqs.(2.5) to other correlation functions.) 
While Tl'^(k 2 ) will have poles corresponding to the pseudo-Goldstone bosons coupling 
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to the currents, Il^ 6 (/c 2 ) is a pole-free function (at least for momenta below the masses 
of the pseudovector resonances). Ilj?(0) is therefore not expected to have a substantial 
dependence on the pseudoscalar masses. The absence of a massless pseudoscalar pole in 
(5.2) requires Uf(0) = 11^(0), so clearly 

i^V^Lo = n ^(°) = n T(0) (5-3) 



Although we do not expect a strong mass dependence of the correlation functions 
(5.1), the individual correlators in the decomposition 



d | 



Anldao^bo-^Wee 



fc=0 



k=0 



d 



+ 2n f 8 b0 M ac -^W S ce 



d 

+ 2nfd a0 M bd — W dS d 
d 

k=o dk 2 55 



fc=0 



(5.4) 



certainly do have an explicit dependence on m. This of course cancels in the sum, although 
the pattern of cancellations is very intricate. 

To illustrate all this, it is instructive to write a simple effective action which 
encodes the information in the zero-momentum chiral Ward identities, and use this to 
derive explicit expressions for the correlation functions in eq.(5.4). 

The zero-momentum Ward identities are 



®acTQ<f>c 



k=0 



2n f S, 



a0 



fc=0 



M ab 



(5.5) 



The simplest effective action compatible with these identities is 

T[QM = Jdx ^Q 2 + 2n/Q*oaV5 + ^"V^-fl/rVs 



(5.6) 



The final term is written in matrix notation. f ab and /i 2 b are matrices, a and f ab are 
constant, and fx 2 b is defined by 



fac l^cd fdb — —M ac & cb 



(5.7) 



(j, 2 is the pseudo-Goldstone boson mass matrix in the OZI limit of QCD, i.e. neglecting the 
coupling to the anomaly Q. 

It is important to realise that the effective action (5.6) is only an approximation, 
where the simplest choice of kinetic terms for the fields <f>^ has been made. This corre- 
sponds to the pole dominance approximation in standard PCAC (see Appendix B). In this 
approximation, as we now show, there is strictly no m dependence in (5.1). 

The second derivatives of the effective action are 



QQ 

bp 



( a" 1 2nf%b 



(5.8) 
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The correlation functions are found by inverting this matrix. We find 



where 



and 



Wee = 


-a A" 1 




W eS b = 






W s * = 


2n f a $ ac (j{k 2 


\ — 1 


Wga gb = 


"$ac A J <S> db 




A = 


l-An)a(f(k 2 - 


"">')» 


A 




4n 2 a I 



-1 



(5.9) 

(5.10) 
(5.11) 



with I = 5 a oSbo- Notice the highly non-trivial fx 2 dependence in all these correlators. (3) 

A simple calculation now confirms that (recalling D a = 2nfS a0 Q + M ac 4>%) the zero- 
momentum Ward identity 

= -M ac $ cb (5.12) 



Wqa qb 
D D 



k=0 



is satisfied. A slightly trickier calculation also shows that the M and /j 2 terms cancel 
completely in the first moment and we find 



d 



Woa 



dtf rvs ° s ° 



k=0 



— facf( 



cb 



(5.13) 



In fact, this is required by the key identity (3.9), which shows 



d 
dk 2 



Wqa qb 
D D 



k=0 



d r 



k=0 



db 



(5.14) 



The first moment of the correlation function W S a S b is therefore given by the decay con- 
stants in the simple effective action (5.6). 

Notice however that / & is only constrained to be a constant by the pole-dominance 
approximation, not by the chiral Ward identities. In general, the / a & appearing in eq.(5.8) 
could be functions of momentum, i.e. f a b(k 2 ). In that case, (5.13) is replaced by 



d 

—rr^Wqa qb 

dk 2 d d 



k=0 



d 

dk 2 V v ' " J k=o 

/ 2 (o)-/(o)^r(o)-/'(o)^ 2 /(o) 



(5.15) 



fc=0 



can have a dependence on the quark masses, but 



We see therefore that W S a^ s ^ 
only proportional to the derivative /'(0). If, as is consistent with the success of PCAC, 



( 3 ) The expressions (5.9) are (in a different notation) identical to those derived in ref. [28] from an effective 
action or diagrammatic resummation including effects beyond leading order in the 1/N C expansion. 
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the pole- free dynamical function f(k 2 ) is slowly- varying in the small- momentum region, 
then this dependence is relatively weak. 



The conclusion of all this is that FF T = W S a S b 



has only a weak residual quark 

fc=0 

owly varying dynamical function. 



mass dependence proportional to the derivative of a s 
The strong, explicit dependence on the quark masses cancels amongst the four individual 
correlation functions in eq.(5.4). 



dk^ w s a D s b D 



The same pattern of cancellations is observed if we now write spectral sum rules for 
using simply the pole-dominance approximation for the propagators given 

/. o 

in eq.(5.9). At first sight, we might therefore expect this pattern to be reproduced in the 
full QCD spectral sum rules described in the next section. In fact that is not so. Of 
course the sum rules go beyond the pole-dominance approximation, but in addition, for 
Green functions such as those needed here, which satisfy dispersion relations requiring 
subtractions, even the sign of the corrections to pole-dominance is not simply determined 
on general grounds. Nevertheless, as we shall see, the full spectral sum rules do confirm 
the picture outlined here by showing numerically the relative insensitivity of -^W S a^ S b 
to the strange quark mass. 



D D 
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6. QCD Spectral Sum Rules and Current Correlation Functions 



First, we derive the QCD spectral sum rules (for a review, see e.g. ref.[29]) for the 
flavour-singlet current correlation function 



(wf 1 1 d ' x etkx{olT d " J ^ {x) dVJ » M |0) (fU) 



Recall the anomalous current conservation equation (c.f. eq.(2.2)) is 

d"J% t = D° = 2m q q l5 q + 2n f Q 

q=u,d,s 



(6.2) 



From its analyticity properties and asymptotic behaviour, the correlator obeys the 
subtracted dispersion relations 



1 

¥ 



dt 



1 



1 



and 



1 

¥ 



Mk z ) - ^5(0) - feVs(O) 



o t (t — k 2 — ie) 7T 

00 dt 1 1 



(6.3) 



lm^ 6 (t) (6.4) 



/o t 2 (t-k 2 - ie) 7T 

From the leading large K 2 = —k 2 > behaviour of the correlator, which is K 4 log(K 2 / /j, 2 ), 



one can deduce that the corresponding derivatives 

f 
(dK 2 ) 



(dK 2 ) 2 \K 2 ) y dK 2 \K*J 



dK 2 \K 4 

are superconvergent and thus obey the homogeneous RGE 

d „ d s d 



dt 



da. 



dxi 



(F;G){t, a s ,Xi) = 



(6.5) 



(6.6) 



Here, Xi = rrii/fj, is the ratio of the renormalised quark mass with the MS'-scheme sub- 
traction scale fj,, and t = L/2 where 



L = log(K 2 /^ 2 ) 



(6.7) 



P, 7 and 7 m are respectively the QCD /3-function, the anomalous dimension for J^ 5 and 
the mass anomalous dimension. The anomalous dimension 7 is 0(a 2 ), viz. 



(6.8) 



and does not contribute at the order we consider here. The coefficients of (3 and 7 m are 



PM = $^&(~) 7m(« s ) = $^7*(— ) 

i=l i=l 



(6.9) 
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where, for three flavours [29], 

Pi = -9/2 (3 2 = -8 

7i = 2 72 = 91/12 (6.10) 

The expression for the running coupling to two-loop accuracy can be parametrised 
as [29]: 

a 8 jjl) 
TV 

where 



a s (p) = ^± = a?) [l - a<°> | log log g + O(aJ) j (6.11) 

= - - 1 . ... (6.12) 
-ft log(/VA) 

and /?j are the coefficients of the j3 function given above. We shall use, for three flavours, 

A = (375 ± 75) MeV (6.13) 

from r decay [30] and LEP[31] data. 

The expression for the running quark mass in terms of the invariant mass m t to 
two-loop accuracy is (4) [32,29]: 

= m, (-/3ia s (^))- 7l//3l [l + |(^-^)a s (^)+0(^)] (6.14) 

where 7$ are the coefficients of the quark-mass anomalous dimension given above. In this 
analysis, we shall retain only the strange quark mass and neglect m u and ma- We use 

m s (l GeV) ~ (197 ± 29) MeV (6.15) 

from the e + e~ — > hadrons[33] and r decay [34] data, and the correlated values of the 
invariant mass m s and A. We shall also use[29]: 

(ss) ~ (0.6 ~ 0.8) (uu) (uu) = -(0.238 GeV) 3 (6.16) 

Now, applying the inverse Laplace operator [35] 

(K 2 ) n d n 

lim (-l) n 7^ TTTTT^T- (6-17) 

K2,n^oo; n/K2=r V ' (fl - 1) ! (8K 2 ) n V ' 

to the dispersion relations (6.3) and (6.4) gives the sum rules ( 5 ) 

r" 3 £(^) + V> 5 (0) = / -j e~ tT -Im^t) (6.18) 

Jo t n 

*£{S) ~ V 5 (0)t + ^(0) = / % e-^-lmfait) (6.19) 

Jo t n 



( 4 ) The truncation of the series at this order is necessary for self-consistency as the quark-quark correlator 
will be used to 0(a s ). 

( 5 ) The unsubtracted sum rule, which is independent of ?/>5(0) and i[>' 5 (0), is more sensitive to the higher 
meson states, and thus is more appropriate for studying gluonium parameters [36]. 
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We use the usual duality ansatz for parametrising the spectral function: 

^1mif> 5 (t) ~ 2m*,f*,6(t - m\) + 6{t-t c ) "QCD continuum" (6.20) 
where f v > is the RG non-invariant rj' 'decay constant' (see refs. [11,37]) normalised as 

(0\d»J% ) \ V ')=2n f V2 f v ,m 2 v , (6.21) 

Then, after transferring the QCD continuum contribution to the lhs of eqs.(6.18) and 
(6.19), we obtain the sum rule: 



V>5(0) 



T 



-C(Gc) + 



1 



rm 2 , 



C{F C ) +^ 5 (0)r 



1 + 



rm 2 v , 



(6.22) 



where the index c in JF C and Q c indicates that the QCD continuum effect has been trans- 
ferred into the QCD expression of the correlators ( 6 ). 

Analogously to JF and Q, their Laplace transforms also obey an homogeneous RGE 
[38], where the resummation of the log-terms can be done by subtracting at r = I /fx 2 and 
introducing the running coupling and masses. 

At this point, it is convenient to decompose the full correlation function (6.1) as follows 



where 



Mk 2 ) - ^(fc 2 )+2(^)^(^ 2 ) + (^) 2 ^(fc 2 ) 

i; gg (k 2 ) = X (k 2 ) = ijd 4 x e lkx (0\T Q(x) Q(0)|0) 

i; qg (k 2 ) = ijd 4 x e ikx (0\T Q(x) s 75 s(0) |0> 

d A x e ikx (0\T s^ 5 s(x) s 75 s(0) |0) 



(6.23) 



ip qq (k 



(6.24) 



Here, we have used the expression (6.2) for the divergence of the current and neglected the 
u and d quark masses. 

For the next stage in developing the sum rules, we need the perturbative QCD ex- 
pressions for these correlators. These have been calculated in the literature. First, the 
perturbative QCD expression for the gluon-gluon correlator ip gg (k 2 ) has been reported in 
[12] and reads: 

iM* 2 ) = + (6-25) 

where 



1+ (^)(; w+ » + . +w+ iA^) 



(6.26) 



( 6 ) In the following, we shall omit this index for convenience of notation. 
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The new corrected coefficient of the perturbative 0(a^) term comes from the erratum 
in [39]. Notice that we have computed the m s -dependent contribution of O{o i s ) coming 
from the 2-loop Feynman diagram with a quark loop inserted on one gluon propagator. 
This is a new calculation that has not previously been published in the literature. (7) 

We have also included the correction due to a tachyonic gluon mass A, where[40]: 

A 2 ~ -(0.43 ± 0.09) GeV 2 (6.27) 

in order to take into account the summation of the perturbative series (as a phenomeno- 
logical alternative to renormalons). A full discussion of the motivation for including this 
term and its phenomenology is given in ref.[40]. 

The non-perturbative contributions come from [41] and [42]. Throughout the analysis, 
we shall use the values of the gluon condensates [43]: 

(a s G 2 ) = (0.07 ± 0.01) GeV 4 (6.28) 

and [41] 

( g s G 3 ) = (1.5 ± 0.5) GeV 2 (a s G 2 ) (6.29) 

A similar value of (a s G 2 ) has been obtained recently from lattice calculations [44]. 

The QCD expression for the quark-gluon correlator ifj qg (k 2 ) has been evaluated in 
[45]. In the MS'-scheme, it reads: 

iM* 2 ) = + (6-30) 

where 

and 7s = 0.5772 ... is the Euler constant. 

Finally, the QCD quark-quark correlator tfj qq (k 2 ) is known to order O(a^) for the per- 
turbative term [46]. Including the condensates of dimension 6 [35,29], it can be expressed 
as: 

^ qq (K 2 ^-k 2 )=K 2 f:^l (6.32) 

d=0 

In the MS'-scheme, the perturbative expression of the renormalised two-point function can 
be written as: 

*>(a- 2 ) = ^£(v) ! C' 33 ) 

i=0 j=0 
( 7 ) We thank Alexei Pivoravov for checking this result. 



' QQ 



lb NP = 
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where c^- are constant terms from the evaluation of the QCD Feynman diagrams. To 
0(a 2 ), they read [47,45]: 



coo 

ClO 



-2 

-131/12 



-17645/24 + (353 - 8w/)C(3) + (511/18)n/ + (3/4)C(4) - 50C(5) 



1 

6 

c i = 1 
en = 17/3 
1 

Cl2 = --jllCQl 

C2i = 10801/144 - (39/2)C(3) -n f [65/24 - (2/3)C(3) 
1 r 

C22 = 7 



cn(-/3i + 271) + 272C01 



1 



C23 = 77^71 (-A + 27i)c i 



where ((n) are the Riemann zeta functions with 

((2) = 4- C(3) = 1.202. 



(6.34) 



(6.35) 



Given the approximations and accuracy to which we are working, we only keep the following 
chiral symmetry breaking and non-perturbative condensate contributions: 



4>2 
V>4 



3 2 



2L + 1 (^) (-3L 2 + 2L - 3 + 6C(3)) 

1 



8tt 2 



:(3-2L) -m s (ss) + — (a s G 



A 112 

V>6 = -m 8 {gs<T' iV -±G° v s) + —7rpa s (ss) 2 



The mixed condensate can be parametrised as: 



(6.36) 



(gsa^'-^G» u s) = M*{ss 



(6.37) 



where the value of M 2 = (0.8±0.1) GeV 2 comes from the baryon [48,49] and B*-B [50] sum 
rules and p = 2 ~ 3 [43,51] indicates the deviation from the vacuum saturation estimate 
of the four-quark operators. 

To complete the sum rules, we need the Laplace tranforms of T and Q. These can 
be obtained from the renormalised QCD expressions above with the help of the generic 
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formulae: 
£ 

£ 

£ 

£ 



where 



1 


n 
T n 


K 2n 


r(n) 


L 




K 2n _ 


- r(n) 


L 2 Z 




K 2n _ 


7(n) 


L 3 Z 




K 2n _ 


~ r(n) 



L 2 + 2iJj(n)L T + V> 2 (n) - ip'(n) 



in 
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It is then convenient to write the Laplace transforms in the sum rule (6.22) in the notation 

C(F;Q) = £(r ga ;Q ag ) + 2(^)£(F qg ;Q qg ) + {^f C{T qq: G qq ) (6.40) 

where the indices gg, qg and qq correspond respectively to the gluon-gluon, quark-gluon 
and quark-quark correlators in eq.(6.24). 



6.1 The sum rule for %'(0) in the chiral limit 

In the chiral limit, ips(0) reduces to the purely gluonic correlator x'(0) evaluated for 
zero quark mass. This estimate has already been done in our earlier paper [NSV2]. In this 
case, the chiral Ward identities require -05 (0) = 0. The Laplace transforms of the gluonic 
correlator read [12]: 
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and 
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(6.42) 



Substituting into the Laplace sum rule (6.22), we find the result shown in Fig. 1. This 
gives a plot of ^/%'(0) versus r for the optimal value of t c = 6 GeV 2 . We observe good 
stability in the range of r from 0.2 to 0.5 GeV -2 . At the stability points, we find 



V^)L.=o = = (26.4 ± 4.1) MeV (6.43) 

The central value is a little higher than the Laplace sum rule value of (22.3 ± 4.8) MeV 
obtained in ref.[12] due to the change[39] in the O(a^) perturbative coefficient in eq.(6.26). 
The different sources of errors are summarised in Table 1. 

Since the validity of the spectral sum rules for calculating %'(()) has been criticised in 
the literature by Ioffe[16,17] (see refs. [16-20] and our detailed rebuttal in Appendix D), we 
should emphasise some features of this derivation. 

First, one should notice that the optimization of the sum rule is obtained at the scale 
t -1 = (2 — 5) GeV 2 , which is relatively high compared to the scale of ordinary mesons of 
the order of m 2 ~ 0.6 GeV 2 . This result is in agreement with the expectation [41,40] that 
the scale of the U(l) channel (gluonium) is relatively high compared with the flavour non- 
singlet (meson) scale. At the optimisation scale we therefore expect (and find) that higher 
dimension condensates (including 'instanton' effects[18,52]) are strongly suppressed. This 
is contrary to the claims in ref.[19], where r _1 is taken at the too low value of 1 GeV 2 . 

Second, the apparently large perturbative radiative corrections in the expressions for 
the two-point correlators tend to cancel in the sum rule (6.22), explaining the almost 
equal value of x'(0) obtained at leading order and the one including the perturbative a s 
corrections. This is reassuring in view of the unknown higher order radiative corrections. 
In order to study the convergence of the perturbative series, we have estimated the a 2 
corrections a la BNP [30] assuming that the coefficient of a s grows geometrically, which, 
numerically, is about 430. This effect remains a small correction to the lowest order 
estimate. 

We have also studied the effect of a 1/k 2 correction due to the summation of the 
perturbative series, which we have parametrised here (see ref.[40]) through a phenomeno- 
logical tachyonic gluon mass A 2 . We see that the presence of this term tends to shift the 
optimisation scale to smaller r values, but affects only slightly the value of x'(0)- 

The decay constant f v > defined in eq.(6.21) can be estimated in the same way using 
just the first (once-subtracted) sum rule in eq.(6.18). We find the result (see Fig. 2) 

/„, = (24.4 ± 3.6) MeV (6.44) 

where the different sources of error are again given in Table 1. 
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One can also derive a Finite Energy Sum Rule (FESR) for x'(0). This can simply 
be found by taking the small r limit of the Laplace transform sum rule. Including the 
radiative corrections, and neglecting the small A 2 corrections, we obtain [12] 



X'(0) 



rf>=<"-(s) 2 ^[ i+ (£) (83 - 4/3i) 



(6.45) 



which confirms the consistency of the set of values of the parameters x'(0), f v / and t c 
obtained from the Laplace sum rules. 



6.2 The sum rule for ^(0) with massive quarks 

In the case of massive quarks, ip5(0) is non-zero. However, its value is known from 
the exact chiral Ward identities (2.9) or (2.11): 
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where we again set m u = ma = 0. As before, we then have two unknown physical 
quantities, viz. limps (t) and ^(O), to be determined from the two subtracted sum rules 
(6.18) and (6.19). (The unsubtracted sum rule[29], independent of ^(O), is more sensitive 
to the higher meson and gluonium mass than the subtracted sum rules, but is in any case 
not needed in this analysis.) 

We now need the Laplace transforms of the quark-gluon and quark-quark correlators 
(see eqs.(6.31) and (6.32)). First, for the quark-gluon correlator, we have 
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The Laplace transform of the perturbative part of the quark-quark correlator including 
the 0(a s ) corrections has a non-trivial analytic expression in terms of the Riemann ( 
functions. It is more convenient to express it in a numerical form: 
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The Laplace transforms of the m 2 and A 2 corrections are 

3 
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and the Laplace transforms of its non-perturbative part are 
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(6.51) 



Finally, collecting all these expressions in the combination (6.40), we deduce the 
Laplace sum rule for the complete correlation function ip' 5 (0) with m s ^ 0. Our result 
is shown in Fig. 1. We find 

^/Vs(0) = (33.5 ± 3.9) MeV (6.52) 

for the range of stability in r = (0.2 — 0.4) GeV -2 . The result obtained with A 2 = is 
very similar, though the stability in r occurs at a slightly larger value of r ~ 0.6 GeV -2 . 

Similarly, we find the value for the decay constant f v > for non-zero strange quark mass 
(see Fig. 2): 

f v > = (27 A ±3.7) MeV (6.53) 

Comparing these results with those obtained above in the chiral limit, we find that 
the effect of the SU(3) breaking quark mass is to increase the values of f v > and -^/ip' 5 (0) by 
approx. 10% and 20% respectively. This is a reasonable conclusion. The SU (3) breaking 
effects are not negligible, but they are of the order expected in relatively smooth quantities 
such as decay constants which, as explained in section 5, are expected to be only weakly 
dependent on the quark masses. 

Certainly we find no evidence of the huge SU(S) breakings advertised in ref.[19], which 
were taken as an indication of the failure of the spectral sum rule method in the flavour 
singlet channel. The main reason for our different conclusion is that the r stability region 
in our calculation is found to be much lower than that used in ref.[19], and so the SU(3) 
breaking terms of (9(m 2 r) are much smaller implying a much better convergence of the 
corresponding OPE. A more detailed comparision of our work with ref.[19] is given in 
appendix D. 
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6.3 Flavour non-singlet correlation functions 



For the full unified Goldberger-Treiman relations, we also need the corresponding 
results for the octet current. These results are immediately obtained from the formulae 
presented above, with the obvious changes of m v for etc. In the derivation, we find it 
convenient to use the normalisations: 

(0|<9<V» = ^| V2/„m;> ( 6 - 54 ) 
^f{k 2 ) = i J d 4 x e ik ' x (0|T D\x) D 8 (0)|0> (6.55) 



and: 



^ 8 (k 2 ) = — ifd A x e lk - x (0|T D°(x) D 8 (0)\0) (6.56) 
2n f J 

where D a is normalised as in eq.(2.10). With these normalisation factors, the quark corre- 
lator i/j ss occurs with the same coefficient in the sum rules for ips, if; 88 and ip® 8 . The octet 
decay constant is related to the conventionally normalised f v by 

f v = 2V6f v (6.57) 

so that /„ = f n = 93.3 MeV for exact SU(3). We find: 

^5 88 (0) = (43.8 ± 5.0) MeV (6.58) 



and 



which corresponds to 



f v = (30.0 ± 3.4) MeV (6.59) 



/„ = (147 ± 17) MeV (6.60) 



The sources of the errors are again tabulated in Table 1. 

In order to quantify the systematic errors of the approach, we have re-estimated the 
value of f n using the same inputs, approximations and methods (for different estimates 
of f n from sum rules, see e.g. [29]) as used above for f v , by using obvious changes of the 
parameters (quark masses, meson mass, continuum threshold). In this way, one obtains: 

/tt = (107 ± 12) MeV (6.61) 

Taking into account the slight deviation of the central value from the experimental number, 
we can consider as a final result: 

f v /U = 1.37 ±0.16 (6.62) 
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where we expect that the error quoted here has been over-estimated. This ratio is in line 
with our expectations, since phenomenologically fx — 1-2/^, and we would expect 577(3) 
breaking to be stronger for the r\ than the K [23,24]. 

As a by-product, we have checked that a systematic rescaling of the value of the 
decay constants f v and f v > will affect similarly the value of the slope of both the non- 
singlet V ; 5 88 (0) and singlet ip' 5 (0) correlators, such that the ratios of correlators and decay 
constants which we use are not affected by this change. 

Together with good t c and r stability, these results therefore confirm the general 
reliability of the spectral sum rules in this channel. 
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Table 1: Error estimates in MeV for the different observables. 
The sources of errors are in [GeV] d where d is the corresponding dimension. 

As already mentioned in sections 3 and 4, the calculation of the off-diagonal correlator 
ip'® 8 is much more delicate, since the contributions of both the r\ and rj' to the spectral 
function must be taken into account. These tend to cancel because of the relative signs of 
the decay constants after flavour mixing, and the sum rule prediction for if;'® 8 should then 
be relatively small. This is confirmed by preliminary estimates. However, a complete cal- 
culation including the effects of flavour mixing for both the off-diagonal i(j'® 8 and diagonal 
V>5 and V'5 88 correlators requires significant further analysis and is beyond the scope of the 
present work. 

We conclude that the QCD spectral sum rule method is indeed giving reliable results 
for the decay constants and susceptibilities in both the flavour singlet and non-singlet 
channels. Although the pattern of cancellations of quark mass effects observed in the 
effective lagrangian, or pole dominance, analysis in section 5 is not manifest in the more 
precise sum rule method, the essential observation that the slope at k = of the correlation 
functions W S a S b is relatively insensitive to the quark masses is confirmed by the numerical 
results found nere. 
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7. Quantitative Analysis of the Unified GT Relation 



In this section, we discuss the implications of these numerical results for the unified 
Goldberger-Treiman relations and the 'proton spin' suppression. 

Collecting the results of the previous section and re-converting to the normalisations 
in sections 2-5, we have found 

(FF T ) = lim -jL W S a Sb = ( (201 ± 23)2 MeV 2 (7.1) 

v >ab fc=0 ^2 d d y (152 ± i7) 2 y v ' 

in the a, b = 0, 8 sector, within an approximation where we have kept only the rj', or 77, in 
the spectral functions for the correlators (0|T D° D |0), or (0|T D 8 D 8 \0), respectively. 

Within this approximation, neglecting 577(3) flavour mixing in the unified GT formula, 
the generalisation beyond the chiral limit of the suppression formula for the singlet axial 
charge is (see eq.(4.10)): 



a 







1 Foo 1 v 7 ^) 



VQ F 88 V^vW^M 



= 0.55 ± 0.02 (7.2) 



This should be compared with the corresponding result in the chiral limit, using our new 
determination (6.43) of x'(0) and using the sum rule estimate (6.61) for 

^8 = 7^ = 0.60 ± 0.12 (7.3) 

The relatively small error in eq.(7.2) is due to the cancellation of the systematic errors 
in the ratio (see Table 1), which was not taken into account in eq.(7.3). Running these 
results from the scale r _1 ~ 3 GeV 2 to the SMC scale of Q 2 = 10 GeV 2 , and substituting 
a 8 = 3F — D = 0.58 ± 0.03, we find 

a°(Q 2 = 10GeV 2 ) = 0.31 ± 0.02 

TP (Q 2 = 10GeV 2 ) = 0.141 ± 0.005 (7.4) 

compared with 

a°(Q 2 = lOGeV 2 ) = 0.33 ± 0.05 

T P (Q 2 = lOGeV 2 ) = 0.144 ± 0.009 (7.5) 

in the chiral limit. 

We therefore find very good agreement between the final prediction for the singlet 
axial charge in the presence of quark masses and in the chiral limit. This confirms our 
theoretical expectation that a is relatively insensitive to the quark masses. Moreover, 
our new prediction for the 'proton spin' suppression remains, notwithstanding the large 
errors on the experimental data, in good agreement with the experimental results quoted 
in section 4. 
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8. Conclusions 



In this paper, we have extended our previous analysis of the 'proton spin' problem in 
the chiral limit by considering the effects of chiral 577(3) symmetry breaking and flavour 
mixing due to the quark masses. 

The formal basis of our analysis is the derivation of the new, unified Goldberger- 
Treiman relations: 

G°a = Fab F v b NN (8.1) 

where F is determined from 

F ac Fl = l^J[LiJdx e ik * (0|T d»j; b {x) <9"4 5 (0)|0> (8.2) 

Apart from the mixing between the flavour singlet and octet channels induced by the non- 
vanishing strange quark mass, the most significant change is the generalisation of the slope 
of the topological susceptibility x'(0) to the equivalent correlation function involving the 
total divergence of the singlet axial current, viz. ^p-(0|T <9 M J° 5 d"J® 5 \0)\ k _ Q . This is the 
quantity which displays the smoothest approach to the chiral limit, with the explicit quark 
mass dependence present in the individual correlators (section 5) cancelling in the sum. 
This observation may have important implications for attempts to calculate the topological 
susceptibility using lattice methods[53], where it is notoriously difficult to approach the 
chiral limit too closely. 

We emphasise again that the unified GT relations are new results and are exact in 
QCD. The familiar PCAC forms (in the flavour non-singlet channels) are obtained by 
approximating the 1PI vertices by the corresponding low-energy meson-nucleon coupling 
constants and by approximating the slopes of the current correlation functions by decay 
constants. Away from the chiral limit, both approximations assume pole-dominance of the 
matrix elements and correlation functions by the pseudo-Goldstone bosons. 

The relevant correlation functions were then evaluated using QCD spectral sum rules. 
In the singlet as well as octet channel, very good stability was obtained with respect to the 
parameter r. The optimum value of r _1 = 2.5 — 5 GeV 2 is in line with general arguments 
for the appropriate scale in the gluon-rich Ua(1) channel. Specific criticisms[16-20] of the 
applicability of spectral sum rules to the Ua(X) channel were shown to be incorrect. 

Our final numerical results agree with our earlier findings[12]. The suppression in the 
flavour singlet axial charge, driven by the mechanism of topological charge screening, is 
confirmed. The reduction in the predicted value of a°(Q 2 ) compared to its chiral limit 
is slight (less than 10%), in agreement with our general theory, although it should be 
emphasised that our results still neglect flavour mixing. Translated into a prediction for 
the first moment of g\, under the assumption that the RG-invariant vertices in eq.(8.1) 
are well approximated by their OZI values, our final result: 

a°(Q 2 = lOGeV 2 ) = 0.31 ± 0.02 

TP(Q 2 = lOGeV 2 ) = 0.141 ± 0.005 (8.3) 

remains in good agreement with experiment, and confirms our proposal that the 'proton 
spin' suppression is a target-independent effect due to the screening of topological charge 
by the QCD vacuum. 
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Appendix A: Partial Legendre transforms 



We derive here the relations between the Green functions and vertices, defined as 
functional derivatives of the generating functionals W and V respectively, used in the 
derivation of the GT relations [11]. 

For generality, we derive these results for a set of fields <E> a and 'currents' J r , with 
sources S a and V r , where the partial Legendre transform is made wrt the sources S a only. 
That is, 

w[v,s] = r[v,$] + s a $ a (Ai) 

We adopt a compact notation where any Lorentz indices are implicit and a spacetime 
integration is assumed in the sum over repeated indices. As in the text, functional differ- 
entiation is indicated by subscripts. Thus, e.g., Ws denotes ^ at fixed V, etc. 
By definition, 

$ a = W s « (A2) 

while 

rv = s a rv = w V r (A3) 

From 

^ = = 6& Wvrsb + 6& Wscsb 



= -r$a $ c Wgcgb (A4) 

we recover the usual result that the 2-point vertex matrix is just the inverse of the propa- 
gator matrix, but in this case restricted to the S a , <E> a sector. 
Similarly, from 

5 , SV* Tjr 5S C Tjr 

5V r 5V r v s 5V r s s 

= Wyrgb — Yyrc&c Wgcgb (^--5) 

we find 

Tyr^b = W V r Sc W~} sb (A6) 
which can therefore be identified as the matrix element (0| J r |$ b ). 



Taking this further, 

o= 6 



WgagcT^c^b j 



5V r 

= Wyrgagc TfycQb + Wgagc (^yr$c$b + Wyrgd T $d <J, c <J, b ^ (A.7) 

from which we find 

Tyrfycfyb + Wyrgd T^dfyc^b = Wgagc Wyrgcgd Wgjgf, (^-8) 
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which is identified as (<£> a | J r |<£> 6 ). 

As a final example, we derive the crucial identity (3.11) used in our derivation of the 
GT relations: 

5 5V f SS d 

Tyry3 = w^ Wys = w^ Wytys + W^ Wsdya 

= Wyrys + Tyr^d Wgdys 

= Wyrys - Wyrgc W~} sd Wgdys (A.9) 

the last line following from eq.(A.6). 



Appendix B: Current algebra, Dashen's formula and the GT relation 

As an illustration of how standard current algebra (PCAC) relations arise from our 
formalism, we give here a short derivation of the Dashen formula for the masses of the 
pseudo-Goldstone bosons. 

The starting point is the identification (see eq.(A.6)) 

T v ^ b = (0\J^\ V b )=tk,r b (k 2 ) (B.l) 

where the on-shell function f ab (m 2 ) is the decay constant matrix for the pseudo-Goldstone 
bosons rf . In terms of the normalised fields rf = B a b4>\, the Ward identity (2.14) is written 

as 

iMV. b „6 + §a C B T cd I\ v - M ac B- b x = (B.2) 

Now take ^f2-| fc _ of this equation. Using the normalisation condition (3.7), we find 
immediately 

/ a6 (0) = $ ac S c T 6 (5.3) 

The assumption that f a i is only a slowly-varying function of k 2 then allows (B.3) to be 
identified with the decay constant matrix. This is (as explained in detail in ref.[ll]) the 
standard PCAC approximation, equivalent to pole dominance of correlation functions by 
pseudo-Goldstone bosons. 

In the same approximation, we can write 

T varjb = k 2 5 ab - (m 2 ) ab (BA) 

Neglecting the mixing with Q, this would produce a pole in the propagator matrix with the 
corresponding mass matrix obtained from the zero-momentum limit of the Ward identity 
(B.2). In matrix notation, this means that 

= -$B T m 2 - MB' 1 (5.5) 
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and so 

$B T m 2 v B$ T = -MB- 1 B<S> T (B.6) 
Using the identification (B.3) of the decay constant, we therefore have (recall <E> T = $), 

fac (m 2 v ) cd fj b = -M ac $ cb (B.7) 

This is Dashen's formula for the pseudo-Goldstone bosons, which for the rif = 2 chiral 
symmetry breaking pattern SU(2)l x SU '(2) _r) SU '(2)y is simply 

flml = -{m u + m d ){qq) (B.8) 

Mixing with the glueball field Q in general shifts the masses of the physical pseudo- 
Goldstone bosons. However, in the simple case where $ a o = 2(<p a ) is zero for o ^ 0, this 
is only relevant in the singlet sector. There, the formula (B.7) refers to the unphysical ry , 
not to the rj'. The corresponding formula for the physical rj' including mixing with the 
glueball field Q, and involving a different identification of the rj' decay constant, is derived 
(in the chiral limit) in ref.[ll]. 

We can also relate the generalised GT relations in the text to the standard formula 
using pole dominance. For rif = 2, eq.(3.18) reduces to 

m N g A = F g nNN (B.9) 

where gA = 2G\ * n our notation. Evaluating F from eq.(3.20) using pole dominance, we 
find 

recovering the standard GT formula and confirming the interpretation of F in eq.(3.18) 
as a decay constant matrix. Notice that since this derivation assumes pole dominance, it 
is an approximation. Unlike the new relation (3.18), the standard GT relation becomes 
exact only in the chiral limit. 



Appendix C: Renormalisation Group 

The renormalisation group equations (RGEs) for the various quantities arising in the 
derivation of the unified GT relations can be derived using the same methods developed 
in ref.[ll]. In this appendix, we summarise the most important identities. 

The starting point is the definition of the renormalised composite operators [54]. For 
QCD, with non-zero quark masses, we have (denoting bare operators with the label '5' 
and renormalised operators with no label for simplicity) 

t0 _ y jOB _ ja^OB 

Q = Q B - 1 ±- f (l-Z)d»J°i 

$1 = z^f r = z^ aB (c.i) 
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where is the inverse of the mass renormalisation, = Z~ 1 . The anomalous dimensions 
associated with Z and Z^ are denoted 7 and 7^ respectively. These definitions ensure that 
the combinations <9 M J° 5 — 2nfQ and m q [q-y^q\ occurring in the Ua(X) anomaly equation 
(e.g. eq.(6.2)) are RG invariant. 

The fundamental RGE for the generating functional W is therefore (in the notation 
of section 2, where suffices on W denote functional differentiation): 

VW = 7(V M ° 5 - W v o h + 70 (s£W s « + S a W s «) +... (C.2) 

where V = ( u-M- + [3-M — 7 m Y] n m Q -J?— ) and a spacetime integration is assumed 

V afI 09 q om q j v,0,S B ,S 

in the sum over repeated indices. The notation + . . . refers to the additional terms which 
are required to produce the contact term contributions to the RGEs for n-point Green 
functions of composite operators. These are discussed fully in ref.[ll], but will be omitted 
here for simplicity. They vanish at zero momentum. 

The RGEs for Green functions are found simply by differentiating eq.(C.2) wrt the 
sources. Simplifying the results using the chiral Ward identities (2.5), we find a complete 
set of RGEs for the 2-point functions. These are: 

VWyO yO = 2'jWyO yO + . . . VW V V b = 7WV0 Vb + . . . VW V a V b = + . . . 

VW v o e = 2 1 W v o. + 7 J-M &W v o sb + ... 

MO MO 2fl f ^ 5 5 

VW V \e = iW V \e + i y-^—M ob W v o S b + . . . 

mj 2nf Vs b 

VWy0 5 Sb = ( 7 + ^)W V 5S b + ... VWy^sb = l^Wya^b + ... 

VW ee = 2-yWee + 2 1 ^—M oh W dSb + ... 

In f j 

VW eS b = (7 + ^)W 9S b + 7-^- {M Qc W S o S b + $ ob ) + • • • 

VW S a S b = 2 l4> W S a S b + ... (C.3) 

It is straightforward to check the self-consistency of these RGEs with the Ward identities 
(2.5) and (2.11). The pattern of cancellations which ensures this is nevertheless quite 
intricate. 

Next, we need the RGE for the generating functional of the 1PI vertices. This follows 
immediately from its definition in eq.(2.12) and the RGE (C.2) for W: 

»r = 7 (v° 5 - J-r Q d„)r v o s - 7 « («r* ; +ri>) + . . . (c.4) 

where V = („£ + ^ - 7m £, m,^) | 
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The RGEs for the 1PI vertices are found by differentiation and, using the Ward 
identities (2.14) to simplify the results, we find in particular: 

VT QNN = -^T QNN + 7- $06 (^Q^Tqnn + TqqT ^j^j + . .. 



2n f 



(C.5) 



Here, V = V + 7<^(</ ,a )^r- As explained in ref.[ll], this is identical to the RG operator 
V defined above (acting on W) when the sources are set to zero and the fields to their 
VEVs. With some calculation, the consistency of the GT formulae (3.4) and (3.5) can now 
be shown using the RGEs (C.3) and (C.5). Again, a very intricate pattern of cancellations 
occurs to ensure this. 

The RGEs for the 2-point vertices occurring in the Ward identities (2.14) are easily 
found by differentiating eq.(C.4). The most important is 

vv n4>l = -274X^6 + (r^Q ($ 0c r^^ - M ob ) + a <-> b ) + . . . (C.6) 

This allows us to deduce the RGE for the matrix B ab which relates the fields to 
the canonically normalised r\ a fields by rf = B ab 4>\. Recall the definition (3.9): 

V 4>Z4>1 1 fc=n = B ac -77^ r 77<V I k=n Bdb = B ac B cb (C.7) 



From eq.(C.6) and the zero-momentum limit of the Ward identities (2.14), we deduce 

VB ah = - l4> B ab + 7 S ac $ c0 $ - 6 1 (C.8) 

The RGE for F ab now follows immediately from its definition F = <&B T and the RGE 
= 7^,$. It is simply 

VF ab = 7 d a0 F ob (C.9) 

that is, 

VF ob = 7^06 for all b 

VF ab = a ^ (C7. 10) 

The final step in proving RG consistency of the unified GT formulae is to show that 
the vertices T v a NN (at k = 0) are RG invariant. Eq.(C.lO) then ensures the required RGE 
for the axial charges, viz. 

VG a A = 1 8 aQ G a A (C7.ll) 

To check this explicitly, notice that eq.(C.5) for T^nn simplifies at k = 0. The contact 
terms vanish and using the zero-momentum Ward identities (see eq.(2.14)) we find 

VT^nn I k=Q = -<y^T^a NN I fc=Q + <y$~Q$ ob r^b NN | fc=Q (C.12) 
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Since r^jvjv 
we confirm 



k=0 




and comparing eq.(C.12) with the RGE (C.8) for B, 




k=0 



= 



for all a 



(C.13) 



Appendix D: Comparison with the literature 

In a number of papers and lectures (see e.g. refs. [16,17]), Ioffe has criticised our earlier 
work on the 'proton spin', expressing his view that our formal theory is "not justifiable" 
and claiming that the spectral sum rule technique we use is not valid for calculating the 
topological susceptibility or decay constants in the flavour singlet channel. 

These criticisms are made explicit in the recent paper [17]. Describing our derivation 
of the Goldberger-Treiman relation (in the chiral limit) for G° A , Ioffe states that 

"the matrix element (p\Q\p) was saturated by contribution of two operators Q and singlet 
pseudoscalar operator <p§ - and the result was obtained by orthogonalisation of the cor- 
responding matrix." This simply reflects a failure to understand our theoretical method. 
As we have repeatedly emphasised, no approximation is involved in the decomposition of 
the matrix element into composite propagators and 1PI vertices, as e.g. in eq.(3.5) here. 
If a different basis of operators is chosen, the definition of the vertices changes too - they 
become 1PI with respect to the new basis. The basis of operators chosen for the decompo- 
sition is in no sense required to be 'complete' - the procedure is not the familiar quantum 
mechanical one of inserting a complete set of states. The matrix element is not being 
saturated with a restricted number of operators chosen from some complete set. The only 
approximation comes in our subsequent conjecture that the vertices, defined specifically as 
we have defined them, obey the OZI rule. The motivations and a posteriori justifications 
for this conjecture are explained carefully and at length in our papers. 

Still referring to our work[12], Ioffe continues, "the calculation of x'(0) by QCD 
sum rules is not correct, because as shown in ref.[19] by considering the same problem 
with account of higher order terms of OPE than it was done in [12], the OPE breaks 
down at the scales characteristic of this problem." This refers to ref.[19], where Ioffe and 
Khodzhamiryan suggested that the extent of SU(3) breaking which occurs in the Laplace 
sum rule approach is unrealistically large, and concluded that QCD spectral sum rules 
were unreliable in the Ua(X) channel. In the rest of this appendix, we show explicitly how 
their calculation is related to ours and point out a number of problems (errors in the QCD 
expressions, inconsistencies of the input and stability parameters, etc.) in their approach 
which are responsible for this false conclusion. 

We consider the current correlation function 



The notation follows our eq.(5.2). Comparing with the notation of ref.[19], the II ^ are 



where J,,l 
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the same up to a minus sign, whereas our Ht is a linear combination of the form factors 
defined there. Only 11^ plays a role in what follows. 

Taking the divergences, and using the chiral Ward identity (2.9), we have 

FrnJ(fc) = k 2 U° L q (k 2 ) = i J d 4 x e lk - x (0\T D°(x) D q (0)\0) + 4m q (qq) (D.2) 



where 



D° = 2n f Q+ 2m qQl5q 

q=u,d,s 

D q = 2Q + 2m q q lb q (D.3) 
The relation with the correlators studied in the text is therefore 

J q=u,d,s 

with the other SU(3) combinations giving t(j® 3 and if}® 8 . Taking the derivative wrt k 2 and 
evaluating at k = 0, we find 

J' q=u,d,s 

lf^{k 2 ) obeys an unsubtracted dispersion relation 

The corresponding Laplace sum rule is simply 

5 £(jr n ) = / dte- tT -ImU° L q (t) (D.7) 
Jo n 



t~ 3 > 



We can also write a once-subtracted sum rule, which enables us to calculate II° 9 (0): 

r°° rli 1 

r- 2 £(^n)+n°/(0) = / ^e-^-Imn^(t) (D.8) 

Jo t n 



Here, 



Taking the SU(S) singlet combination of these sum rules for q = u,d,s clearly gives just 
the Laplace sum rules (6.18) and (6.19) analysed in the text. 
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In ref.[19], Ioffe and Khodzhamiryan define the 'decay constants' for individual quark 
flavours as follows: 

(o|jyv> = *M?' (£-io) 

and by saturating the r.h.s. of eq.(D.7) with the rj', obtain sum rules for f v f gl> for q = u,d 
and q = s. Taking the ratio, they find an unrealistic 577(3) breaking characterised by 
gL/gp ~ 2.5 and conclude that the Laplace sum rule for II ^ with the flavour singlet axial 
current does not work. 

However, as we have shown in the text, a correct implementation of the sum rule 
method does indeed work, showing good stability over a range of appropriate r and t c 
values and giving results for the decay constants in both flavour singlet and non-singlet 
channels which show only the expected level of 577(3) breaking. In the rest of this section, 
we shall therefore quote the full formula for the sum rule (D.7) for II^ 9 , correcting some 
mistakes and omissions in ref.[19], then briefly indicate some of the problems with their 
calculation. 

From the results given in section 6, we can immediately read off the required sum 

rule: 

2n /V ^/„<m;,e-™S' + ... = ^r" 2 (^) 2 (l - (1 + t c r)e~^) (1 + 6%) 

3 (OL s 



+ i(v)l^< a « G2 >( 1 + ^ + (v)^ G3 > 



«.^"X(l-e- tcT )(l + 
2(1 + 35 qs ) (^)M$Tm a {ss) - 5 qs 4m s (ss) 



(7J.11) 



where (see also [35]) 



u gg 



( ^^ +/3l(1 _ 7£)+6 ( 3m! + ^) T 

9tt\ 



XNP 

6 99 =~ 



TV 



J qq 



(D.12) 



The dots on the l.h.s. denote the presence of any further intermediate states (in particular 
the rj) which might be required. The terms appearing with 5 qs are present only when q is 
chosen to be the s quark, and as usual we have assumed m u = ma = 0. 

We now comment on some of the differences between the paper [19] and our work: 
• It is a priori dangerous to use the sum rule for with the individual quark flavours 
while only keeping the rj' as an intermediate state in the sum rule. This would imply 
that only the rj' and not the rj would contribute in the mixed SU(3) combination n° 8 . It 
follows that the quantity g^,/g™l d used in ref.[19] is not the most reliable measure of the 
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strength of SU(S) breaking: while a value close to 1 would indicate self-consistently that 
SU(3) symmetry is very accurate and the 77 — 77' mixing angle is very small, a large value 
would simply show that the initial hypothesis of 77' saturation for the individual U ^ is not 
self-consistent. (Ref.[19] refers to a further calculation taking 77 — 77' mixing into account, 
but the details are not given.) 

• The formula (D.ll) differs in a number of important ways from that used in ref.[19]. 
In (D.ll), the coefficient of the term 4m s (ss) coming from the zero-momentum Ward 
identities has been corrected compared to ref.[19]. We have included the leading-order 
radiative corrections, which in this example have a significant effect in inducing stability 
in t and t c . Notice that the radiative corrections have a far greater effect on this single 
sum rule than in the combined sum rule used in the text to calculate ^5(0) , where they 
cancel between the C{F) and C(Q) terms. We have also included the contribution of 0{m 2 s ) 
arising from the correlator denoted by ifj qq in the text, which was omitted in ref.[19] but 
is of the same order as the other principal terms for the values of r used in [19]. 

• In our analysis of the sum rules, we find good stability at values of t c ~ 6 GeV 2 and r in 
the range (0.2 — 0.4) GeV -2 . In contrast, ref.[19] uses too low a value of t c ~ 2.5 GeV 2 and 
without finding stability takes an ad hoc r value of about 1 GeV -2 . This is far too high 
and goes against the expectation [41,40] that the optimal scale r should be comparatively 
low in the gluon-rich flavour singlet channels. For this value of r, the convergence of the 
OPE for the terms of (9(m 2 r) is also poor. 

We can therefore conclude that the analysis of the sum rule for Ti^k 2 ) in ref.[19] 
(equivalent to our once-subtracted sum rule for ips(k 2 )) is unreliable and that their con- 
clusion that the spectral sum rule method does not work in the flavour singlet channel is 
mistaken. 

By taking the once-subtracted sum rule (D.8) for T\°®{k 2 ) in combination with (D.7), 
we can write a Laplace sum rule for n° (0), the quantity required for the generalised 
Goldberger-Treiman relations (and therefore the 'proton spin' effect). This sum rule is of 
course identical to that for ^(O) analysed in the text. This is not discussed in ref.[19]. 
Alternatively, Ioffe and Khodzhamiryan use a finite energy sum rule (FESR), analogous 
to the one described in our eq.(6.45) but omitting the radiative corrections. However, the 
FESR should be used with great care due to its strong t c dependence, and the fact that 
the result comes from a difference of large numbers. Indeed, it is unsurprising that by 
substituting their set of values of f v >, g v ', and t c , which are inconsistent with the values 
we have obtained from the Laplace sum rule based on r and t c stabilities, they obtain an 
unreliable result. 

Finally, in more recent papers[20], Ioffe and Oganesian extend the work of ref.[19] by 
combining it with a sum rule for G° A itself, derived using a composite 3-quark operator with 
the quantum numbers of the proton. However, one can immediately notice that the choice 
of the interpolating nucleon currents used by the authors is far away from the optimised 
choice analysed in refs. [48,49]. Irrespective of the validity of this method, it is evident that 
the results of ref. [20] are as unreliable as those discussed above, and for the same reasons, 
since the results of ref. [19] are used as input parameters in the new sum rule. 
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Fig. 1 The dependence of the correlation functions x'(0), y/i/j' 5 (0) and ^/ ip'® 8 (0) (in GeV) 
on the Laplace sum rule parameter r (in GeV -2 ). 
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